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We study the corrections to the quenched limit of ^ > QCD. We use an improved reweighting procedure. 



Introduction In Fig. |l| we sketch tentative fea- 
tures of the QCD phase diagram in the k, /x, T 
space. Since a general description is missing, one 
may try to concentrate on special situations, such 
as the "quenched" limit (see below), repre- 
sented by the thick vertical dashed line at infinite 
\x in the k = plane. For a more physical sit- 
uation we consider a vertical plane at small but 
non-zero k, then, e.g., a possible transition would 
appear at /ic(K, T) < oo. Here we present a first 
attempt to describe phenomena away from the 
quenched limit. 
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Figure 1. Tentative phase diagram. 

Hopping parameter expansion The QCD 

grand canonical partition function is (we specify 
for one fiavour of Wilson quarks and neglect con- 
stant factors; IJ are links, T lattice translations): 

Z(/3,«,/i) - 1 [i?C/]e-^«(^-^^»Z;^(Ai,/i,{C/}), 
Z;^(/t,/.,{C/}) = DctW^(/t,M, {[/}), (1) 
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U,T, + r^,T*U*) 



(e ^ r+4 C/4 T4 + e F_4 TT W 



with r±p = 1±7^, 2k 



1 

m+4 ' 



(2) 
Here 



m-t-3-t-cosh u 

m is the naive continuum limit bare mass" , m 
the bare mass at /i = 0. The analytic expansion 
in K leads to an expansion in closed loops, the 
links contributing factors k or ne , see (^) : 

ZFiK,fi,{U}) = DctT4^ = cxp(Tr lnT4^) 



exp 



EEE 

i=l {C,} s=l 



n n , 

{C,} 



(3) 



1=1 



(D, C = Dirac, colour). Ci are distinguishable, 
non-exactly-self-repeating closed paths of length 
^, s is the number of times a loop Cci covers Ci, 



fci 



ee 



Y if = ^^Polyakov r—path'\ (4) 



1 otherwise. A ^^Polyakov r-^atK^ closes over the 
lattice in the ±4 direction with winding number 
r and periodic (antiperiodic) b.c. [e = -|-1(— 1)]. 
Quenched limit at /i > In the double limit ||] 
K ^ 0, /i ^ 00 with Ke'* = C fixed, only straight, 
positive Polyakov loops Vg are retained: 



z]?l(C,{;7}) = exp 



Ef:^Trc(7',)^ 



{x} s=l 

eCVs)\ C=(2C) 



(5) 



- \{ Detc(l 

{x} 

(see II , and for further studies). C (or C) 
fixes the direction toward the limit - see Fig. 1. 
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Next order corrections To disentangle the ef- 
fects of K, and /i we go to next order in n: 



{x} s=l 



{ecy 



Tr 



x,r,q,i,t,t' 



r,q 



(6) 



For easy bookkeeping we use the temporal gauge 
Un,4: = 1, except for U(^s,ni=N^)A = V^: free, and 
-P^Lt' = (V5)'-~^t/(^,),(14+i)«t/(^,,),, (7) 
with r > q>0, i ^ ±1,±2,±3, 1 <t <t' < 
(t < t' for g = 0). As an approximation to the 
full theory one may expect to be good in some 
small region around the vertical limit line and by 
varying fi or T we may hope to get information, 
say, about the phases around /ic at large mass. 
Alternatively we may propose as a model by 
itself. At large /i and mass it would be an ap- 
proximation in the above sense at all T. At large 
T it would also be an approximation for any ji, as 
we recover the "Polyakov loop model" there Q]. 
We may thus hope that our model retains some 
features of QCD with full k, h dependence. 
Improving convergence of /j, > simula- 
tions Since at > Zp is complex, simu- 
lations using reweighting methods suffer from a 
"sign problem". In this situation one can try to 
realize the cancellations between contributions in 
a correlated way. This can be done by identify- 
ing appropriate transformations of the configura- 
tions and correspondingly separate the measure 
into symmetric and antisymmetric factors, to be 
used in the MC generation and reweighting, re- 
spectively. In special cases this method can be 
developed to a solution of the sign problem |^ , in 
our case we may hope to achieve some improve- 
ment - see also Q. The sensitive dependence of 
Zp on the temporal Polyakov loops suggests to 
use center rotations. This is especially sugges- 
tive for our model (||-^), but the method can be 
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Figure 2. History of weight (above) and density 
(below), iV^ = 4, /3 = 5.5, n = 0.12, /i = 0.6267. 

applied also to the full problem. We write 
Z = y"[DC/W]e-^^«^'^»Zi.({i7,y}) 

= j[DUDV]S{{U,V})T{{U,V}) (8) 

with center-symmetric V"}) = S{{U, ZV}), 

Z = {l,e2W3 e4W3}^ and calculate 



(O) 



X{{U,ZV}). (9) 



Here are averages taken with the Boltzmann 
factor S. U denote spatial links. We choose 



-Sg 



Y^Zf{{U,ZV}) 



(10) 



Alternatively, in the approximation of retaining 
only the s = r = 1 terms in (^ we can take for S 



(11) 



Note that S already includes k and /i. In Fig. g we 
show the convergence of the weight T and charge 
density averages (|2|) for each Z in (||) and after 
symmetrization. Note the good convergence in 
a regime of large cancellations (C = 0.0407; the 
weight is only ~ 3%). The symmetrization im- 
proves the statistics by a factor 3-10 in the con- 
fining phase. In the deconfining phase, whichever 
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sector is chosen by S the method automatically 
captures the right contribution. For the analysis 
we use 8'^4 and 8^6 lattices with e = —1, perform- 
ing ^ 10^ heat-bath/over- relaxation/metropolis 
sweeps with @ or (|l]) and r = 1 (first 20000 
sweeps for thermalization). Errors are estimated 
with a jackknife analysis. 
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Figure 3. Above: ns/T^ vs C, various k. Below: 
ub/T^ at fixed C and ms/lO vs k^. 

Results We calculate the baryon number density 
ub and the "mobility" niB (we introduce this as a 
normalized measure for the effect of the term): 

T-^ no + ni 



no 



= 2C(^TrPs), 



(12) 
(13) 

(14) 



as well as ((X^jTrP^)^), the spatial and tempo- 
ral plaquettes and the topological susceptibility. 
The effect of the correction can be observed 
in Fig. H. As expected, the mobility depends pri- 
marily on K, not on /x, and the correction to the 
density depends linearly on k"^ and it can be sig- 
nificant. The most important point, however, is 
that we can now probe the k and /i directions in- 
dependently. The question of phase transitions is 



investigated in Fig. ^, both along the T and along 
the ^ direction. We observe the installation of the 
deconfining regime with increasing T (the tem- 
poral Polyakov loop is non-zero in the confining 
phase due to the fermions). Also, a signal for a 
transition in /i is seen, especially at small temper- 
ature. We cannot say, however, from these pre- 
liminary results whether these are genuine phase 
transitions or mere crossovers. The further devel- 
opment implies a more detailed analysis. 
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Figure 4. Above: nB/T^ and Polyakov loop vs 
13. Below: nB/T^ vs ^. 
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